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ABSTRACT 

Recently, some intriguing results have lead to speculations whether the central density slope - velocity dis- 
persion anisotropy inequality (An & Evans) actually holds at all radii for spherical dynamical systems. We 
extend these studies by providing a complete analysis of the global slope - anisotropy inequality for all spheri- 
cal systems in which the augmented density is a separable function of radius and potential. We prove that these 
systems indeed satisfy the global inequality if their central anisotropy is (3q ^ 1/2. Furthermore, we present 
several systems with (3q > 1/2 for which the inequality does not hold, thus demonstrating that the global 
density slope - anisotropy inequality is not a universal property. This analysis is a significant step towards an 
understanding of the relation for general spherical systems. 
Subject headings: galaxies: kinematics and dynamics - dark matter - methods: analytical 
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1. INTRODUCTION 

Theoretical dynamical models continue to play a key role in stellar dynamics, as understanding their underlying structure helps 
shed light on the properties of numerical and observational stellar systems and dark matter haloes. In this paper, we focus our 
attention on the relation betwee n the density slope y (r) and the velocity anisotropy profile (3(r), which has attracted renewed 
interest lately. As is well known. lAn & Evan s (2006) proved that the central inequality 7n ^ 20q is a necessary co ndition for the 
positivity of the distribution function (DF) of a spherical system. More recently however Ciotti & Morganti ( 2010a b) showed that 
-f(r) 2/3(r) at all radii (hereafter called the Global Density Slope - Anisotropy Rela tion, GDSAI) i s a necessary c ondition for 
positivity of the D F, if /3 ^ 1/2, in the families o f multi-component Osipko v-Merritt dOsipkovll 1 9791; [Merritt 1985), Cuddeford 
Cudde ford 1991), and Cuddef ord-Louis models ([Cud deford & Loufsl ll995l) . as well as for the Plummer models oflDej onghe 
1987), the Hernquist models of Baes & Deionghd (12002ft . and the models we introduced in Baes & van Hesd (120071) (hereafter 



Paper I). Their proof is based on the fact that all these models are characterized by having a separable augmented density (see 
Section|3]l. They also note that currently, no counter-examples of the GDSAI are known, but remark that in the case of Cuddeford 
models with a central anisotropy (3q > 1/2 the GDSAI is only a sufficient condition, so that possible counter-examples could be 
found in this range of values. 

These results pose the question under which conditions the GDSAI holds for all spherical systems. In this paper, we make 
important advancements by providing a complete analysis of the GDSAI for all well-behave d systems with a sepa rable augmented 
density. This group includes all aforementioned models, as well as the hypervirial models of | Evans & Art (|2005). the 7-models of 
iBuvle et all (120071) and the Dehnen-McLaughlin systems discussed in lVan Hese et al.l d2009l) (hereafter lPaperlll) . among others. 
Fir st, we show that the GDS AI holds for all separable systems, if (3q ^ 1/2, by proving an equivalent criterion formulated 
by lCiotti & Morgantil d2010b"t) . In this manner, we extend their previous results. Our analysis also reveals some very peculiar 
properties of separable systems. Furthermore, we show that counter-examples of the GDSAI do exist for separable systems with 
(3q > 1/2, in other words, we demonstrate that the GDSAI is not a universal property. However, the velocity distributions of 
these models are extreme, and all counter-examples are very likely dynamically unstable. 

First, we outline in Section [2] the general concepts of spherical dynamical models. In Section [3] we describe the augmented 
density framework. In Section|4] we give our analysis of the GDSAI for separable systems: we prove the inequality for models 
with /?o ^ 1 /2, and we present three counter-examples with /?o > 1/2. Finally, we discuss our results in Section|5] 

2. SPHERICAL DYNAMICAL MODELS 

The dynamical structure of a spherical gravitational equilibrium system, governed by a positive potential ip(r), is completely 
determined by the non-negative phase-space distribution function (DF) F(r, v). For spherical systems, this DF is a function 
F(E, L) of the isolating integrals, the binding energy E and the angular momentum L: 

E = W)-\vl-\v%, (1) 
L = rv T , (2) 

with 

v t = \Jvl+vl, (3) 

the transverse velocity. From the DF, the velocity moments 

M2n, 2m (r) = 2ttM J j F(E, L) v 2n v 2 T m+1 dv r dv T - (4) 
can be obtained, with M the total mass of the system. In particular, the density and the second-order moments are 

p(r) = fioo(r), P&r( r ) = /i 2 o(r), pc4(r) = (5) 
and cr|,(r) = 2a 2 (r). The density slope and the velocity anisotropy profile are defined as 

7(r) = -^(r), (6) 



Spherical dynamical models satisfy the Jeans equation 



dlnr 



ftr) =l-irl (7) 



which can be written as 
with 

Evidently, it follows that 



dpa 2 2p(r) 2 di/f 

( r ) + — — P°r (r) - p(r) — (r) , (8) 

a 2 r (r) ( 7 (r) - 2/3(r) + «(r)) - v 2 (r), (9) 

, . dine;: , . 2 , , dip . . 

K(r) = -— —(r), v 2 c (r) = -r-^(r). (10) 
d in r dr 

7 (r) - 2/3(r) + K (r) > 0, Vr. (11) 



ICiotti & Morgantil (l2010 a.b) showed that several systems (see Introduction) satisfy a stronger condition, the GDS AI 

7(r) - 2/3(r) ^ 0, Vr, (12) 

and they pose the question whether this condition holds for all spheri cal systems. Natur ally, the inequality is valid outside the 
radius r 2 where j(r 2 ) = 2. It is also valid at r — 0, as was proven bv lAn & E vans (2006). In this paper, we will investigate the 
GDSAI for a particular class of systems, namely those with a separable augmented density. 

3. THE AUGMENTED DENSITY CONCEPT 

A spherical dynamical system can also described by an augmented velocity moment (Deionghe ( 1986); Paper I), which extends 
the moment to an explicit function jj.2n.2mi4>, t) of both the radius and the gravitational potential. An augmented moment is 
equivalent to the DF: the knowledge of one augmented moment determines the entire system. In particular, we will consider the 
augmented density p(ip, r), and its relationship with the DF is given by 

, x f* fW- E ) F(E,rv T ) 2 

p(!l>,r)=2irM AE J> dv 2 . (13) 

Jo Jo V 2(V> - E) - vf 

This integral equation can in principle be inverted to obtain the DF by using Laplace-Mellin transforms, although in practice the 
inversion is only numerically stable for sufficiently smooth systems. The strength of the augmented density framework lies in its 
direct connection to observable quantities like the velocity moments. For instance, the augmented velocity dispersion profiles are 
given by 

v 2 M, r) = -^-^ / W, r) <ty', (14) 



p{ip,r) 
2 



■4, 



tT(4>,r) = —— D r , r 2 p(V',r) <ty', (15) 



p{ip,r) 

where D r i denotes the derivative with respect to r 2 . The observed density and dispersions are then simply recovered from 

p(r) = p(ip(r),r), (16) 
a 2 r (r) = a 2 (i>(r),r), (17) 
a 2 (r) =a 2 T {i>{r),r), (18) 

and the density slope is a ~ ■ , a ~ 

, > r op , , . . r aw , , op , , , 
7(r) - —jrO>(r),r) - --f r ^^ r ,r . (19) 
p or p or o%p 

As remarked in the Introduction, Ciotti & Morganti have examined the GDSAI in several systems with a separable augmented 
density, i.e. systems of the form 

~p{iJ,r)=f{4>)g{r), < V < V>o, (20) 
with V^o = 4>(0). For such models, the dispersion profiles read 

^={ 1+ l^)m)J? m ^ ,22) 

Note that the radial velocity dispersion is now only a function of i/j. The velocity anisotropy profile of these systems has the 

simple form , , , 

= --— — (r . (23) 
2 d In r 

As we demonstrated in Paper I and Paper II, this property provides a very elegant way to construct dynamical models with a 
given potential, density and velocity anisotropy. Indeed, separable systems are completely determined by tj}{r), p(r) and (3(r), 
since g(r) is defined by Eq. d23l and, by inverting t/j(r), the function f(ip) follows from 

fW - P ~P§y (24) 

However, one still needs to verify whether the corresponding DF is non-negative everywhere. Eq. ( Il9l now reduces to 

d In g , , d In ib , . d In f . , . , 
7(r) = ~^(r) - -^(r) —±(i>(r)), (25) 
d In r d In r d In yj 

so that we obtain . r f( I ( \\ 

4(^)) = Z ^(7«-2/3(r)). (26) 
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In other words, as remarked by Ciotti & Morganti, the GDSAI 

7(r) > 2/3(r), Vr > 0, (27) 
is for separable systems equivalent to the statement 

^y>0, VO^V^o- (28) 
dtp 

The question thus becomes whether this inequality is valid for all separable systems. In the following section, we will prove that 
this is indeed the case, if /3q ^ 1/2. 

4. ANALYSIS OF THE GDSAI FOR SEPARABLE SYSTEMS 
Following the reasoning of lAn & Evansl (120061) . we assume that any well-behaved DF can be written in the form 

F(E,L) = L~ 2/3 ° (F Q (E)+F 1 (E,L)) , (29) 

with 

F 1 {E,0)=0, V0^E^tp Q . (30) 

The function L~ 2 P° Fq(E) in this Ansatz can be understood as the leading term of a Laurent series expansion in L at L = 
0. Towards the center r —> 0, the DF is dominated by this term, which has the form of a system with constant anisotropy. 
Consequently, the central anisotropy of the entire model indeed corresponds with /3q. Since the DF has to be non-negative 
everywhere, it follows immediately that Fq (E) ^ WE is a necessary condition to obtain a physically meaningful DF. 
If we consider separable systems, the corresponding augmented density then has the form 

pty, r) = fty) r~ 2Pa (1 + 9l {r)) , with 5l (0) = 0. (31) 

2 

Using u 2 = 2 (^~e) ' me re l a ti° n between the augmented density and the DF dTJb can be written as 

p(ip,r) = 2n2 1 / 2 - 0o r- 2f} °M t " = d,, 2 j - £)V2-A ( F (E) + F x (e, ru^2(ip - E)) 

Jo VI- u 2 Jo V v 

In separable systems, it follows that 

m = ~wr- { } 

Since the left-hand side of this equation is independent of the radius r, the right-hand side does not depend on r either. The 
equality is therefore valid for all values r; in particular, we can take the limit of r towards the center, 

f(tp) = lim P^t-LLL = li m r 2fio p(iP, r). (34) 

r— ¥0 g\T) ' — >0 

This property is the key element to prove the GDSAI when /3 ^ 1/2: using (1301 and (|34l , it follows from Eq. ( 1321 that 

f{tP) = {2-Kf'H-^M r yf~^\ - Efl 2 -'^ F (E)dE. (35) 

1 [6/1 - po) Jo 

Remarkably, the function f(tp) thus only depends on Fq (E) and (3 . In other words, for separable systems the function i*\ (E, L) 
has no influence on the GDSAI. Concrete exampl es of this behavior are fur nished in the systems considered by Ciotti & Morganti. 
For instance, the equivalent function B(tpj) in lCiotti & Morganti (2010b) for generalized Cuddeford systems does not depend 
on the anisotropy radius r a (see their Eq. (13)). 
The value of $ splits our further analysis into three cases: (3 < 1/2, (3 — 1/2, and f3 > 1/2. 

4.1. Proof for /3 < 1/2 

If Po < 1/2, the derivative of f(tp) becomes 



dE. (32) 



= (2nf/ 2 2^M r ^ 1 " /3 °\ 
dip r(3/2-/3 ) 



lim (tp - E) 1 ' 2 -^ F (E) 



+ (\-fa)J (^-E)- 1 / 2 -^F (E)dE 



(36) 



Let us examine the first term inside the brackets: if 



lim (tp - E) 1/2 - p ° F (E) > 0, (37) 

then 

lim (tp- E)- 1/2 - l3a F a {E) ~ lim (tp - E)- a with a ^ 1, (38) 
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so that 



(tp - Ey 1/2 ~P° F {E)dE = +00. 



(39) 



In other words, if the limit is nonzero, then the integral in the second term becomes infinite. The limit can therefore be omitted, 
so that the equation is simplified to 



*L m = fo^a/ao-* T (l ~ ft) f* F (E) 



dE > 0, 



(40) 



and recalling Eq. (126*1 ). the GDSAI is proven. The above relation can be generalized further: if n = [3/2 — /3oJ and a 
3/2 — J3q — n are the integer floor and fractional part of 3/2 — /3q, then 



so the inequalities 



^W = (2vr) 3 / 2 2-^M y-ft>) 
d^ fcV ^ V ' r(3/2 - ft - fc) y 



< fc sc n, 



fc WO > 0, 



(41) 



(42) 



are nece ssary conditions to obtain a separable system with a non-negative DF. This extends the results obtained bv lCiotti & Mor ganti 
(2010a) for multi-component Cuddeford models. 



When /3 Q = 1/2, Eq. OB reduces to 



The derivative is then simply 



4.2. Proof for f3 a = 1/2 

/(VO = 2ir 2 M [ F (E)dE. 
Jo 

^j)=2tt 2 MF q (iJ;)^0, 



(43) 
(44) 



so evidently, the GDSAI is again a necessary condition for a physical dynamical model. 



4.3. Counter-examples for /3 > 1/2 

The proof is not applicable to systems with /3q > 1/2. Indeed, the derivative has the same form as Eq. Pol l, but now the two 
terms inside the brackets are respectively +00 and — 00 when Fq(E) > 0, so their sum is undetermined. However, we can rewrite 
Eq. 03 using integration by parts as 



^3/2-0 o £ Q ( ) + / ty, _ £;)3/2-ft, F^(E)dE 



where Fq(E) denotes the derivative of F (E). After differentiation, we then obtain 



(45) 



(46) 



Thus, separable systems with a monotonically increasing Fq(E) (i. e. Fn(E) ^ Vffl, satis fy the GDSAI. Again, this is an 
extension of the results for generalized Cuddeford systems found bv lCiotti & Morgantil (1201 Obi) . 

Yet, the GDSAI is no longer a necessary condition for a physical model, which raises the question whether systems can be 
found for which the global inequality does not hold. To this aim, we consider the potential-density pair 



ip(r) 
p(r) 



GM tc 



1 + 

3M 
"87 



with corresponding density slope 



7(r) 



r 3 / 2 (l + Vr) 
3/2 + 7/2^ 



(47) 
(48) 

(49) 



1 + vV 

which is part of the family of Veltmann models or a- models (Veltmann ( 1979j);|Zhiil[l996)), and was discussed bv lMoore et al.l 
(119981) . If M to t = M, then the system is also self-consistent. For this pair, we construct physical DFs that generate four-parameter 
anisotropy profiles of the form 

M s ft + fto (r/r*) 2S 



1 + (r/r a ) 
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with < S < 1, so that 

p(^r)=/W(-) [l + ^l , (51) 

with 



,. N -2/3o ( r 26\> 



P<5 = 7 • ( 52 ) 

Again, our systems have separable augmented densities. For every anisotropy profile, the function f(ijj) follows from p(r) — 
p(ifj(r),r), and t he DF c an be found by inverting Eq. ( fT3l ). Instead of performing these calculations directly, we adopt the 
technique used in lPaper l5 we first generate a family of components of the form 



~ Pi ^,r)=poi[±) | 1 — ] [-] + . (53) 



where pi, qi and Si are three parameters, and poi are normalization constants. The corresponding DFs are (see lPaper"H : iPaperll) 

POi ST^f 1M / & \ I & 



-/3 +fc<5 



y (h\ r(i+ P<+J - a< ) / _L^_ | (54) 

^ V k ' r (i - a, + fej) r ( ft + Js . t -i + ^o-M) V 2r ^/ 



fori 2 < 2r.fi, and 



Pi+j Si -3/2 



X 



y- ni+Pi+jsi) / js_ ] (55) 

& v fe / r(i- J 9 0O -*j)r(p i +ja,-| + i 8 0O + fc*) V 2r * 2£ 7 



for i 2 > 2r 2 i. With different values of the parameters, a library of base functions is thus created, from which a linear 
combination is built that fits the given density p(r) at various radii. This is achieved by minimizing the quantity 



m=l \ i=l / 



Xat = ^— > , -7— r i /">»' ' ■>."'/'- ( '•»<' I • i56) 

using a quadratic programming algorithm (Deionghe 1989). The details of this procedure can be found in Paper II. In particular, 
we created models with N = 12 components, fitting 25 density data points extracted from Eq. d48l l. 

With this technique, we obtain several dynamical models with non-negative DFs that violate the GDSAI; three of them are 
shown in Fig.Q] All three share the anisotropy parameters /3q = 0.75, (3^ = 1 and r a = 0.02, but have different values for 5: 0.3, 
0.6 and 1 .0 respectively; note that the latter is a Cuddeford-type model. For the model with S = 0.3, we find that j(r) < 2/3(r) for 
radii in the interval ]0, 0.021], with a minimum around r = 0.0057 (note that the center is a local maximum, for which 70 = 2/3o). 
In the model with 5 — 0.6, the 7 — f3 relation reaches a local maximum around r = 0.0028, and the GDSAI does not hold in the 
interval [0.019, 0.061], with a minimum around r = 0.036. Finally the largest 7 — f3 fluctuations occur in the Cuddeford model 
{5 = 1), with a local maximum around r = 0.0054, and a GDSAI violation within [0.019, 0.100], with a minimum for r = 0.044. 

Evidently, we require rather extreme parameter values to obtain these (modest) violations, while maintaining non-negative 
DFs. The central anisotropy /3q has to be high, and the profile f3(r) has to increase very rapidly. It is therefore safe to assume that 
the self-consistent variants of these models are dynamically unstable. This can be seen from the standard criterion for radial-orbit 
instability: 2T r /Tx = 2(v r ) / (vt) = 5.45, 8.26 and 10.42 for the three models, which is much higher than the ~ 2 threshold 
for similar models (see iMerrittl (11999b for an overview). Further evidence of dynamical instability is given by the radial velocity 
distributions 

F Vr {r) = 2irM / F(E,L)v T dv T - (57) 







As shown in th e bottom row of Fig. \T\ these profil es have two or three peaks at small radii. These are indications of Henon 
instabilities (see IMerrittl (11999b : [Barnes et al] dl986l) ). In theory, if the systems are instead not self-consistent but embedded in a 
massive dark matter halo, they might withstand these instabilities; however, one can safely argue that such equilibrium systems 
are too extreme to arise in structure formation. 
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4.4. The inverse relation 

Finally, we rema rk that the function F (E) can be derived from f(ip) by means of an Abel-related inversion (ICuddefordll991t 
lAn & Evansll2006l) . which holds for all values of ffa < 1, 

F(F] = / r E d " + V d^ j_ £7 A 

01 ; (27r) 3 / 2 A/r(l - a)r(l - /3„) io d^™ +1 " ^ Q d^ n / ' 



where again n = [3/2 — /3 J and a = 3/2 — /3 — nare the integer floor and fractional part of 3/2 — /3 - Thus the additional 

rW>o, V0^^4 (59) 



condition ^n+i j 



is sufficient to obtain a non-negative Fq{E). As lCiotti & Moreantild2010al) showed, this also implies that the entire DF F(E, L) 
is non-negative in the case of generalized Cuddeford systems. However, it is not a priori clear whether this property is true for all 
separable systems, since the behavior of Fi (E, L) might still lead to negative values of the DF. Further study is therefore needed 
to determine if Eq. ( |59l is a sufficient condition for the existence of a physical separable model. 

5. DISCUSSION 

In the previous section, we presented a full analysis of the GDSAI for spherical dynamical systems with a separable augmented 
density. As our proof shows, the GDSAI hold if the central velocity anisotropy /3q ^ 1/2. We further demonstrated that systems 
with fto > 1/2, the GDSAI can be broken, as shown by three counter-examples, although these systems are not physically 
realistic. 

Eqs. ( |40b and d44i l. combined with Eq. (|26*T i, also reveal a remarkable property of separable systems: the GDSAI is purely 
determined by (3q and Fq(E). The function Fq(E) can be interpreted in various ways: it can be thought of as the phase-space 
distribution of particles at purely radial orbits, as the phase-space distribution of particles at the center, or as the energy distribution 
of the constant-anisotropy component of the DF. As a surprising consequence, if a separable system has a given potential ip( r ) 
and density p(r), then knowledge of Fq(E) alone is sufficient to construct the complete DF of the system. Indeed, we showed 
that Fq(E) is equivalent with f(tp), and in combination with p(r), the function g(r) = p(r)/f(ip(r)) can also be derived, 
determining the augmented density f(ip) 9(f) and thus F(E, L). 

The next logical step will be to investigate the GDSAI for general, non-separable spherical models. One possible approach 
would be to consider a spherical systems as a linear combination of separable systems. In fact, an analytic p(tp, r) or F(E, L) 
can be written as a double sum of power-law functions, by means of a two-dimensional Laurent series expansion. An alternative 
approach would be to ask the following question: given a spherical dynamical system with a given ip(r) and a DF that generates 
p{r) and /3(r), does there always exist a separable model with a non-negative DF that generates the same density and anisotropy? 
As we mentioned in Section[3] the function g(r) follows directly from /3(r), and in turn this determines f(ip) from g(r) and p(r). 
However, there is no a priori reason that the corresponding DF is also non-negative. If it is, the same GDSAI analysis applies as 
presented in this paper. We currently do not know of any (p(r), /3(r)) pair that can only be generated by non-separable models, 
but more study is required to resolve these questions. We think our analysis of separable systems can be a useful stepping stone 
for further investigations of the GDSAI for spherical dynamical systems. 
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FIG. 1. — Three models for which the GDSAI does not hold: <5 = 0.3 (solid line), 5 = 0.6 (dashed line), and S = 1.0 (dotted line). In the first panel, the 
density data points are also displayed. 



